Astronomy & Astrophysics manuscript no. 7741 


© ESO 2008 


February 1, 2008 





The nature of the fluorescent iron line in V 1486 Ori 
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ABSTRACT 

The fluorescent 6.4 keV iron line provides information on cool material in the vicinity of hard X-ray sources as well as on the characteristics 
of the X-ray sources themselves. First discovered in the X-ray spectra of the flaring Sun, X-ray binaries and active galactic nuclei (AGN), 
the fluorescent line was also observed in a number of stellar X-ray sources. The young stellar object (YSO) V1486 Ori was observed in the 
framework of the Chandra Ultra Deep Project (COUP) as the source COUP 331. We investigate its spectrum, with emphasis on the strength 
and time variability of the fluorescent iron K^, line, derive and analyze the light curve of COUP 331 and proceed with a time-resolved spectral 
analysis of the observation. The light curve of V 1486 Ori shows two major flares, the first one lasting for ^ 20 ks with a peak X-ray luminosity 
of 2.6- 10^^ erg/s (dereddened in the 1-10 keV band) and the second one - only partially observed - for > 60 ks with an average X-ray luminosity 
of 2.4 ■ 10^' erg/s (dereddened). The spectrum of the first flare is very well described by an absorbed thermal model at high temperature, with 
a pronounced 6.7 keV iron line complex, but without any fluorescent K„ line. The X-ray spectrum of the second flare is characterized by even 
higher temperatures (> 10 keV) without any detectable 6.7 keV Fe XXV feature, but with a very strong fluorescent iron K„ line appearing 
predominantly in the 20 ks rise phase of the flare. Preliminary model calculations indicate that photoionization is unlikely to account for the 
entire fluorescent emission during the rise phase. 
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1. Introduction 

The COUP data set provides a nearly continuous 13 day long 
X-ray observation of the Orion nebula star forming region and 
is currently one of the best suited data sets for the exploration 
of the X-ray properties of large and homogeneous samples of 
young stars. Substantial a nalysis of this data set has aheady 
been performed (e.g. by iGetman et al.l l2005h . In particular, 
Tsuiimoto et alj (120051) carried out a systematic search for the 
occurrence of the fluorescent iron Kq, line at 6.4 keV, identify- 
ing 7 COUP sources with detectable fluorescent emission. This 
emission line is particularly interesting since it can provide a 
wealth of information on the emitting source(s). Fluorescent 
6.4 keV photons - hereafter referred to as Fe I K„ photons - re- 
sult from a de-excitation of neutral or lowly ionized iron atoms, 
following an excitation leaving the iron ion with a K-shell va- 
cancy. When the latter is filled by another (mostly L-shell) elec- 
tron, the excess energy can be carried by a fluorescent photon. 
Highly ionized iron can also show this behavior, but the transi- 
tion energy rises because of the reduced screening of the inner 
electrons from the nuclear charge. Several processes lead to the 
ejection of a K-shell electron and, thus, to the creation of the 
Fe I K„ line. During the process of photoionization a K-shell 
electron is ejected from iron after the absorption of a photon 
with an energy exceeding 7.11 keV. Alternatively, high energy 
particles can eject K-shell electrons with the same threshold 
energy. 



The first detection of the fluorescent Fe I Kg line was repo rted 
in the context of solar physics (e.g. ^ Neupert et al. 1967 ). A 
model succeeding in explaining most of the solar Fe I Kq^ emis- 
sion with p hotoioniza tion only was prop osed and evaluated by 
Bail (Il979b . yet, e.g., ^Z arro et al] (Il992 ) argue that photoion- 



ization alone is not sufficient to account for all of the solar 
Fe I Kq, emission. Fluorescent K^, line emission has also been 
reported from other celestial X-ray sources. It is typically en- 
countered in the high-energy emission from AGN and X-ray bi- 
naries, which are characterized by high X-ray luminosities and 
rather hard X-ray spectra, thus providing large ph oton numbers 



above the 7.11 keV threshold for photoionization. iFabian et al 



(11982) 

detect the fluoresc ent iron line i n the spectrum of the 
X-ray binary Cyg-Xl, and lTanaka et al.l d 1995b were the first to 
report the relativistically broadened line in the spectrum of the 
AGN MCG-6-30-15. Recent observations with XMM-Newton 
and Chandra extended the class of Fe I line emitters to a 
number of stellar s ources, such as the class I YSO Elias 2-29 
dFavata et a l. 2005), the ab ove seven COUP sources discussed 
by 'Tsuiimoto et al. (2005"), and the YSO YLW 16A for which 
Imanishi et al. (2001) carried out time-resolved spectroscopy. 

In this letter we present a detailed analysis of the X-ray 
properties of the YSO V 1486 Ori (= COUP 331) with em- 
ph asis on its Fe I Kq. em ission, which was already reported 
by iTsuiimoto et al.l (l2005i) . who measured a normalization of 
2.6(0.1 - 4.6) • 10~^^ ph/(cm^ s) with an equivalent width 



2 



S. Czesla and J.H.H.M. Schmitt: The nature of the fluorescent iron Une in V 1486 Ori 



of 126 eV for the fluorescent line. Unfortunately, the YSO 
V 1486 Ori has so far not received much attention at other en- 
ergy bands. From the 2MASS-all-sky survey we infer its in- 
frared magnitudes as 12.5 mag, 10.6 mag, and 9. 4 mag, in the 
J,H, and bands, respectively. ITsuiimoto et alJ (l2005h show 
that the infrared colors of V 1486 Ori do not match the re- 
gions of reddened dwarfs or giants and argue for the presence 
of a NIR-emitting inner disk. Moreover a rotation pe r iod of 
(6.0 9 + 0.3) d is known for V 1486 Ori ('Herbst et al.' ('200^ 



COUP 331 2-9 keV light curve 



and IStassun et al. (119991) '). Following Ts uiimoto et al.. (.20051) 



we adopt a distance of 450 pc for Orion. 

2. Observations and data analysis 

The COUP data was obtained during 13 consecutive days 
in Jan. 2003 with the Advanced CCD Imaging Spectrometer 
(ACIS) on board the Chandra X-ray observatory. The total ex- 
posure time is » 840 ks, separated into 6 observation segments 
extending over a total time of ^ 1 140 ks. All observations have 
the same aimpoint (a = 5*35'" 17 ' and 6 = -5°23'40") and the 
field of view covers an area of 17'x 17'; the source COUP 33 1 is 



locate d at a = 5"35'"9.2' and 6 = -5°30'58" (cf.'GetmaneLal 



2005h . Our analysis is based on the pre-processed photon data 
provided by the Chandra pipeline; for further analysis we ap- 
plied the CIAO-software in version 3.4. We first screened the 
photon data for events in the 0.3 - 13 keV energy band to 
suppress background. Thereafter, we defined source and back- 
ground regions as a circle with an 8 arcsec radius and an annu- 
lus extending from ^ 14-32 arcsec, both centered on the nomi- 
nal source position. 

2.1. The light curve 

We generated the background-subtracted 2-9 keV band light 
curve of COUP 331 shown in Fig. [T] The 6 observation seg- 
ments can be clearly identified. For most of the COUP observa- 
tions COUP 33 1 was recorded at a count rate of 3.3 • 10"^ cts/s, 
except for two periods. A strong flare (termed "flare I") oc- 
curred x 90 ks after the beginning of the observation. The flare 
rapidly rises to its peak within 4 ks, the decay is again quite 
rapid with an e-folding time of r^fcay = 5.2 ks, but, unfortu- 
nately, is only partially observed. A second, less prominent but 
very significant count rate enhancement (termed "flare 11") oc- 
curred 1080 ks after the start of the COUP campaign. This sec- 
ond flare was also not fully covered. Flare 11 has a much slower 
rise phase lasting for 20 ks, and thereafter, the 2-9 keV 
band count rate remains more or less constant for at least 40 ks 
when the COUP observations were terminated. For the rest of 
the observations COUP 331 was found in quiescence. 

2.2. Spectral analysis 

In the following section we carry out temporally resolved spec- 
tral analysis of the flare I, flare II, and quiescent phases of 
COUP 331; note that all quoted errors refer to 90% confidence 
intervals unless stated otherwise. COUP 331 is located about 
7 arcmin off-axis, hence it does not appear as a perfectly sym- 
metric point-source. However, we regard its image distortion 
as weak and apply the CIAO standard tools for the extraction 
of point-source spectra (i.e. psextract). The spectra were an- 
alyzed using the XSPEC environment in version 1 1 .3. 1 . 
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Fig. 1. Background-subtracted light curve of COUP 331 with 
2 ks binning. Inserted panels show close-up views of the flares 
I and II. 

2.2.1 . The quiescent phase 

The quiescent phase spans the - 1050 ks time interval with 
the exception of the first flare event, and is the longest last- 
ing phase, covering « 90% of the total observation time. Note 
that we disregard aw 20 ks contribution preceding flare II. 
During quiescence COUP 331 is detected with a net count-rate 
of 3.3 • lO"-' cts/s corresponding to an (unabsorbed) X-ray lumi- 
nosity of 1.6 ■ 10^" erg/s. We checked for temporal changes in 
the quiescent spectrum, but did not find evidence for any sig- 
nificant variations. The quiescent X-ray spectrum is well de- 
scribed by an absorbed thin-thermal plasma model (we used 
the VAPEC models) with all metal abundances - apart from 
Fe, which is left as a free parameter - fixed to 0.3 times the 
solar values. In Fig.|2]we show the data and the fit; the best fit 
model parameters are presented in Table [l] where we list the 
derived absorbing hydrogen column-density, the temperature, 
the iron abundance and the fit quality. 

2.3. Flare I 

For an analysis of flare I we considered all data recorded 
in the time interval from 87.5 to 102 ks; the spectrum of 
flare I is shown in Fig. |2l The flare spectrum is very well de- 
scribed by the same model as applied for the quiescent phase 
(cf. Sect. 12.2. Il l with somewhat diff'erent spectral parameters 
listed in Table [T] We find - as expected for a flare - an en- 
hanced temperature of 6.8 keV, leading to a spectrum with a 
clearly detected Fe line feature at 6.7 keV, but with no detected 
excess emission at 6.4 keV. From the spectrum we derive an 
(unabsorbed) peak X-ray luminosity of Lx = 2.6 • 10^^ erg/s 
(in the 1-10 keV band) reached directly after the rise phase, 
and - with an e-folding time of 5.2 ks for the decay phase - a 
total energy output of AE ^ 1.4 • 10^^ erg in the same band. 
With these numbers, flare I observed on COUP 331 is among 
the largest stellar flares ever observed. 

2.4. Flare II 

For the analysis of flare II we considered all data recorded after 
1080 ks (cf. Fig.[TJ. The X-ray spectrum was generated and an- 
alyzed in the same fashion as those of the quiescent and flare I 
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COUP 331 during quiescence and fiare I - Data and model {absorbed VAPEC) 

' "m,iMlUJf#,,,,l,i 




Energy [keVl 

Fig. 2. Spectrum of COUP 331 during flare I (upper curve) and 
the quiescent phase (lower curve), both fitted with an absorbed 
VAPEC model. 

phases. From the spectrum of the entire observed flare II phase 
(see Fig. [3]) we calculate an 1 - 10 keV band X-ray luminosity 
of 2.4 • 10-" erg/s. Considering a duration of at least 60 ks, the 
total energy output in this energy band is AE > 1 .4 • 10^^ erg. 

In the fit process of the flare II data we encountered two 
problems. First, the absence of temperature-sensitive features 
in the spectrum - especially the hot Fe XXV feature at 6.7 keV 
and its Fe XXVI counterpart near 7 keV - prevents us from 
tightly constraining the plasma temperature. We find that any 
choice of the temperature beyond 10 keV leads to acceptable 
fits, and thus decided to apply a value of 13 keV. Second, 
the spectrum shows definite excess emission above the thermal 
continuum at 6.4 keV, which we attribute to the Fe I K„ line. 
We therefore included an additional instrumentally broadened 
Gaussian component centered at 6.4 keV into our model and 
refitted, leaving only the normalization of the Gaussian com- 
ponent free to vary. The thus obtained fit is shown in Fig. [3] 
which shows the obvious presence of the 6.4 keV line and vir- 
tual absence of the 6.7 keV Fe XXV feature. The resulting fit 
parameters for the thermal model are again listed in Table [T] 
and for the Gaussian component we arrive at a normalization 
of 5(3.3 - 6.7) • 10"* ph/(cm^ s), corresponding to an equiva- 
lent width of 685(452-918) eV. As is obvious from Fig. |3 the 
inclusion of a Gaussian component at 6.4 keV into the model 
significantly increases the fit quality. 

Table 1. Fit parameters for the absorbed VAPEC model. 
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Fig. 3. The COUP 331 spectrum during flare II with a close-up 
of the 5-8 keV energy band. Two model fits are indicated. First, 
(solid line) an absorbed thermal model and second (dashed 
line), the same thermal model with additional instrumentally- 
broadened Gaussian component centered at 6.4 keV. 

narrow -band light curve of COUP 331 with a 4 ks binning in 
addition to the 2 ks binned 2-9 keV band light curve; both 
light curves are shown in Fig. |4l As is clear from Fig. |4]the 
2-9 keV band light curve stays more or less constant after the 
rise to peak, while the 6.1 - 6.7 keV narrow-band light curve 
decays quite rapidly. We therefore subdivided the the flare II 
phase data into three consecutive 20 ks time intervals, denoted 
by "rise-phase", "phase 11", and "phase III" as indicated in 
Fig. |4] The X-ray spectra in the 5-8 keV region correspond- 
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As a first result we thus note that the 6.4 keV emission feature 
is present - at best - only during 60 ks during flare II. In or- 
der to further temporally constrain the presence of the 6.4 keV 
iron line we constructed a background-subtracted 6.1 -6.7 keV 



Fig. 4. The background-subtracted narrow -band light curve 
(solid line) of COUP 331 during the last 80 ks of the COUP 
observation. The light curve peaks in the rise-phase of flare II. 
For clarity the 2-9 keV light curve is also shown (dotted line 
- scaled to fit the plot range). 

ing to these time intervals are shown in Fig. |5] In these plots 
the solid line indicates the thermal model described earlier For 
all spectra the same model was applied, but we allowed the nor- 
malizations of the VAPEC and Gaussian component to be fitted 
independently for all time intervals; all other parameters were 
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regarded globally. Figure |5] clearly demonstrates an evolution 
of the Fe I K^, line strength. Interestingly the line normaliza- 
tion decreases as time proceeds, while the underlying thermal 
spectrum remains almost unaffected. We thus interpret the vari- 
ations in the narrow -band light curve shown in Fig.|4]as arising 
from variations in the Fe I Kq, feature's strength. Note that the 
constancy of the thermal spectrum does not imply the same for 
the temperature, since the latter exceeds 10 keV, producing a 
flat spectrum in the "low energy" band we are observing. The 
Kq. feature is strongest during the flare 11 rise phase where our 
fits yield a line flux of 9. 1 (5 .8 - 1 3 .4) 1 0"^ ph/(cm2 s) (an EW of 
* 1400 eV) while in phase 11 and 111 it is reduced to 3.04(0.8- 
7.6) and 0.9(0-4.8)- lO"*" ph^cm^ s), respectively This is also 
reflected by the narrow-band light curve around 6.4 keV (cf. 
Fig. Si. 
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Fig. 5. Time-resolved 5.5 - 8 keV section of the COUP 331 
spectrum (cf. Fig. [Hi. Note the decrease of the Fe 1 K„ line's 
strength. 

3. Interpretation and conclusions 

Our time resolved analysis of the COUP X-ray data of the YSO 
COUP 331 shows that for ^ 90% of the entire observation 
time the object is found in a quiescent state characterized by 
an X-ray temperature of 30 - 40 MK. Two major flares are 
covered by the COUP observation. In the first data segment a 
relatively short duration event with a peak X-ray luminosity of 
2.6 ■ 10^^ erg/s occurred, which clearly shows the 6.7 keV iron 
feature, but no Fe I K^^ emission. 

In the last COUP data segment a long duration flare event 
occurred; starting with a a; 20 ks rise phase it lasted > 60 ks. 
The decay of this flare cannot be constrained because the 
COUP observations were stopped. The total amount of energy 
released in the 1-10 keV energy band is > 1.4 • 10^* erg; the 
X-ray temperature of this flare is extremely large (> 10 keV) 
and can in fact not be well constrained with the Chandra data. 
A clear and highly significant detection of the fluorescent K„ 
line is obtained only during the 20 ks rise time of this flare; at 
later times the line may still be present albeit at a lower level. 
This event on CO UP 331 resembles a huge flare observed on 
Algol, observed by lFavata & Schmittl(ll999h with BeppoSAX, 
which also shows a slow rise phase, followed by a nearly con- 
stant phase - lasting for w 80 ks - in its 1.6 - 10 keV band 



light curve. However, the corresponding 15 - 100 keV band 
ligh t curve clearly shows the characteristics of a decay phase 
(cf.. ^ Favata & SchmittI (1999), Figs.l and 2) . The energy out- 
put of 1.4 ■ 10^^ erg derived bv lFavata & Sch mitt (1999) in the 
0.1-10 keV band compares well to the numbers estimated for 
COUP 331, note however that no Fe 1 K„ line was detected. 

The precise temporal association between Fe 1 Kq, line 
emission and high energy continuum is necessary in order to 
quantitatively assess the origin of the Fe I Ka emission ob- 
served in COUP 331; we find that the actual photon flux in 
the fluorescent Kq, line during the flare II rise ph ase is 35 times 
larger than reported by Tsuiimoto et al. ' ('2005'). This discrep- 
ancy arises because Tsuiimoto et al. (2005) give time-averaged 
values, whereas our values refer to several limited time in- 
tervals. The X-ray temperature during the flare event was ex- 
tremely large; the total luminosity in the Fe I Kq line alone is 
X 2.3 • 10^" erg/sec. While these circumstances appear to fa- 
vor the excitation of the Fe 1 Kq line through photoionization, 
preliminary model calculations challen ge such a sc enario. We 
repeated the calculations carried out by Bail (Il979h applying a 
reflector with cosmic abundances. We used the same branching 
ratios and cross sections for iron, but slightly different cross 
section for absorption by other elements composed according 
to cosmic abundances. Furthermore, we located the source di- 
rectly on the reflector to make it subtend half the sky, thus pro- 
viding an upper limit for the fluorescent iron line flux. The 
illuminating input spectrum was modeled as a thermal X-ray 
spectrum with a temperature of 13 keV, normalized to the mea- 
sured fluxes. Our preliminary calculations show that an already 
optimally assumed fluorescence geometry significantly under- 
predicts the observed Fe I Kq flux, making an interpretation 
of the Fe I K,, line through photoionization difficult to accept. 
However, these calculations will be extensively discussed in a 
larger and systematic context in a forthcoming paper, where we 
will also address its impact on the characteristics of COUP 331 
and other sources with fluorescent K„ emission. 
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